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(~| ^ Non-Gaussian fluctuations of the electrical current can be detected with a Josephson junction 

I , placed on-chip with the noise source. We present preliminary measurements with an NIS junction 

Y\ ' as a noise source, and a Josephson junction in the thermal escape regime as a noise detector. It 

is shown that the Josephson junction detects not only the average noise, which manifests itself 
as an increased effective temperature, but also the noise asymmetry. A theoretical description of 
the thermal escape of a Josephson junction in presence of noise with a non-zero third cumulant 

"j^ ■ is presented, together with numerical simulations when the noise source is a tunnel junction with 

Poisson noise. Comparison between experiment and theory is discussed. 
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1 Introduction 



The fluctuations of the electrical current reveal the charge of the carriers, their correlations as 
^ \ well as their ferniionic nature. Until recently, only the variance of the current fluctuations had 

been measured and predicted. Theoretical progress in the last decade allowed one to calculate the 
full counting statistics of electrons passing through any phase coherent circuit made of conductors 
f^ ■ such as tunnel junctions, quantum point contacts, diffusive wires or chaotic cavitics[TJ[31[31[31[S]. 

For short enough measurement times, the distribution of the number of transmitted charges is 
predicted to differ sizcably from a Gaussian distribution. On the experimental side, few exper- 
f~^ \ iments have gone further than measuring the average noise, i.e. the quadratic average of the 

^D ■ fluctuations. Indeed, observing a departure from the Gaussian distribution of the fluctuations 

requires the measurement of small signals in a short time with a high accuracy. The pioneer- 
ing experiment performed by Reulet et al. used analog microwave techniques to reconstruct 
{V^) , with V the voltage across a current-biased, ~ 50 rj-resistance tunnel junction. This first 
jrt ' measurement of the third order cumulant of the voltage fluctuations (or skewness) pointed out 

the importance of feedback effects associated with the electromagnetic environment of a noise 
source [71 |B]. Bomze et al. [5] performed a direct counting of single tunneling events across a 
tunnel junction, using a high-precision analog-to-digital converter. Another approach consists in 
using an on-chip detector |10|. The statistics of the charge in quantum dots in the classical regime 
of sequential tunneling was probed with a quantum point contact [TTJ [12]. It was also shown 
that the current through a Josephson junction in the Coulomb blockade regime is sensitive to 
non-Gaussian noise |13] . A promising scheme to probe current fluctuations, proposed by Tobiska 
and Nazarov [14j . consists in using a Josephson junction as a threshold detector. The switching of 
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a Josephson junction out of the zero-voltage state provides a sensitive threshold detector because 
the switching rate varies exponentially with the bias current. We present here a preliminary ex- 
periment in which a current-biased Josephson junction is capacitively coupled to a voltage-biased 
tunnel junction which induces current fluctuations in the Josephson junction. A similar experi- 
ment was performed by Timofeev et al. |15j using a junction in a regime involving macroscopic 
quantum tunneling, thermal escape and phase diffusion. The response of the Josephson junction 
to the squewness of current fluctuations was compared with an adiabatic model which uses as an 
input the measured response to the average current and as a fit parameter a bandwidth, of the 
order of the plasma frequency. In our experiment, the junction is kept in the regime of thermal 
escape which allows direct comparison to theory. The switching rate of the Josephson junction is 
measured for opposite values of the average bias current lb through it, thus giving access to the 
asymmetry of the current fluctuations. In a second part, theoretical predictions for the response 
of the Josephson junction to fluctuations with a non-zero third cumulant are presented, under the 
assumption that the correlation time of the fluctuations is shorter than the response time of the 
Josephson junction, and that feedback effects related to the voltage developing on the Josephson 
junction can be neglected. The results of this calculation are then compared with numerical simu- 
lations that integrate the Langevin equation describing the phase dynamics in presence of current 
fluctuations created by uncorrelated tunnel events. Finally, experimental results are confronted 
with theory. 

2 Measuring noise with a Josephson junction 

In a simplistic picture, a Josephson junction biased at an average current If, below its critical 
current Iq can be seen as a threshold detector for current fluctuations: neglecting thermal activa- 
tion or quantum tunneling, the junction remains in its zero-voltage state as long as the current 
fluctuations are smaller than Iq — Ii,. When the current exceeds this limit, the junction switches to 
the dissipative state, and a voltage develops. This voltage subsists until a "reset" is performed. 
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Fig. 1 Simplified scheme of the experiment. The high frequency (/ > 1/RC) fluctuations Sim 
of the current through a voltage-biased tunnel junction pass through a Josephson junction. The 
switching of the Josephson junction to the finite- voltage state depends on the sum of Sim and 
of the pulsed current bias Ib{t). The switching probability of the junction during one pulse can 
be related to the current fluctuations. 



by driving the bias current down to zero, then back to /&. The behavior of the Josephson junction 
is thus similar to that of an electrical circuit breaker. The principle of our experiment is to couple, 
through a capacitor, the current fluctuations 6Im of a voltage-biased noise source to a current- 
biased Josephson junction (see Fig. 1), and to compare the switching rates r+ and r_ at opposite 
values of the bias current ±Ib applied during a time Tp through the Josephson junction. The 
asymmetry in the rates Rr = r+/r_ — 1 reflects the difference in the probabilities that during Tp 
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fluctuations such that 5Im > lo — lb and Sim < —{Iq — h) occur, i.e. the asymmetry in current 
fluctuations whose amplitude exceeds Iq — lb. 

However, a Josephson junction is not a perfect threshold detector: the switching of a Josephson 
junction occurs for bias currents lb smaller than /g, because of thermal activation or quantum 
tunneling. As explained below, thermal activation is the dominant process in our experiment; we 
will thus consider this case only. Any dissipative impedance connected to the Josephson junction 
produces Johnson-Nyquist current noise, which can yield thermal escape and switching. The 
theoretical description of this phenomenon is based on the Resistively Shunted Junction (RSJ) 
model, in which the dynamics of the Josephson junction is similar to that of a particle in a tilted 
washboard potential with friction [16l[T7j. The state of the Josephson junction is characterized by 
the phase difference ip across the junction, which relates to the position of the fictitious particle. 
The particle evolves in the Josephson potential —Ej{cos ip+sip), where Ej = Iq^Pq is the Josephson 
energy, with (/3g = h/2e, and s = Ib/Io the reduced bias current, which determines the average 
slope of the potential. At small bias current, the particle is trapped in a local minimum and 
oscillates at the plasma frequency ujp{s) — ujo{l — s^)^^^, where the frequency ujq ~ (/o/(/3qCj)^/^ 
is set by the parallel capacitance Cj. When the current increases further, the thermal energy can 
become comparable to the barrier height AU{s) w (4^2/3) (PqIo{1 — s)^'"^. The particle can then 
be thermally activated over the barrier and escape from the local minimum at a rate 
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In this expression, the escape temperature T^sc is a measure of the current noise in the circuit 
connected to the Josephson junction, mostly at the plasma frequency. In the setup presented in 
Fig. 1, two noise sources are connected to the Josephson junction: the resistance R in series with 
the tunnel junction, that can be assumed to be at the bath temperature To, and the tunnel junction 
itself, whose current noise is described by a spectral density 2e (Im) coth(ei?t (Im) /2kBTo), with 
Rt the tunnel resistance. In a simple model, the escape temperature in Eq. ^ is therefore given 
by[l8]: 

Tesc = (l/R + l/Rt)-\To/R + e (/,„) /2kB coth{eRt (Im) /^kB^)). (2) 

In this setup, the Josephson junction is thus very sensitive to the second moment of the current 
fluctuations [II HO]. 

The measurements presented here were obtained on a sample with an aluminum Josephson 
junction of area 0.7 /im^, in parallel with a capacitor Cj = 36 pF. For convenience in sample 
fabrication, the noise source under investigation was not a normal (NIN) tunnel junction, but 
a 0.2 /im^ NIS junction made of Al and Cu electrodes separated by a tunnel barrier obtained 
simultaneously with the barrier of the Josephson junction by oxidation of the aluminum layer. In 
order for the electromagnetic environment of the Josephson junction to be entirely determined 
by on-chip elements, large capacitors (Ci « 850 pF and C sa C2 ~ 130 pF) were fabricated on- 
chip between the connecting pads and ground (see Fig. [2]). These capacitors and the junction 
capacitor Cj were made using plasma oxidized aluminum electrodes deposited on the substrate 
prior to patterning the circuit |18] . At frequencies of the order of the plasma frequency, the 
capacitors act as short circuits in parallel with the circuit placed outside the chip. Capacitance 
Ci also allows for a good voltage bias of the junction and reduces the low-frequency cut-off of 
the high-pass filter associated with Rf and C. Resistors Rf — 202 Vt and i?(, = 202 Vt were 
commercial surface mount NiCr resistors connected directly on the sample with silver epoxy. The 
sample holder was thermally anchored to the mixing chamber of a dilution refrigerator with a 
base temperature of 20 mK. All lines connected to the sample were twisted pairs equipped with 
microwave filters, except a 50 fi coaxial line with distributed 50 dB attenuation designed to apply 
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Fig. 2 Full scheme of the experiment. The noise source is a Cu-Al NIS tunnel junction. 
The current pulses through the Josephson junction are obtained by voltage pulses on a coaxial 
line, applied to the series combination of a bias resistor Rt and the Josphson junction. Large 
capacitors C, Ci and C2 allow for good voltage bias and limit the impedance seen by the 
Josephson junction to on-chip elements. 



current pulses to the Josephson junction. The floating bias of the NIS tunnel junction was realized 
with a battery-operated voltage divider, followed by an electrical two-position switch that allowed 
reversals of the sign of the bias (at ~ 1 Hz) to eliminate the effect of drifts. Bias resistors of 10 MO 
were placed on the two wires of the twisted pair connecting the voltage source to the samples. 
The voltage was measured across the series combination of the tunnel junction (noise source) and 
of the Josephson junction (detector), allowing for measurements of one or the other, depending on 
the bias scheme. For example, the I — V characteristics of the NIS junction shown in Fig. [3] was 
measured with /{, = 0, so that the voltage drop across the Josephson junction was zero. From this 
characteristic taken at 20 mK, we deduced the tunnel resistance Rt — 1.57 kl7. In the following, we 
describe measurements with the NIS junction biased above the gap voltage, in a transport regime 
dominated by tunneling of single electrons, like in a usual NIN tunnel junction. 
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Fig. 3 a) Current-voltage characteristics of the NIS junction at a refrigerator temperature of 
To = 20 mK. A fit to theory yields the tunnel resistance of Rt = 1.57 kil and the superconducting 
gap A = 185.5 fieV.h) Current-voltage characteristics of the Josephson junction at two different 
values of the current /,„ in the NIS junction and at a refrigerator temperature of To — 20 mK. 
The black dots were measured with (Im) = 5.3 nA and the grey dots with (Im) = 640 nA. 
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The current-voltage characteristics of the Josephson junction, shown in Fig. [3] for two values 
of the current through the NIS junction, is obtained by sweeping back and forth the voltage 
Vfc, and recording both the current (from the voltage across a lOfi resistor placed at the base 
temperature) and the voltage. A supercurrent branch is observed, and when Vb/Rb exceeds a 
threshold Isw, the junction switches. The current then drops essentially to zero, with the applied 
voltage Vb = Rblsw '^-^ 0.1 mV < 2A/e now entirely across the Josephson junction. Compared to 
a perfect current bias, our bias scheme limits the production of quasiparticles at the Josephson 
junction. However, we observed in the I — V characteristics of the Josephson junction a large 
sub-gap current peak close to V ~ 0.15 mV, leading to increased heating for the largest switching 
current. We therefore had to reduce the duty cycle of the applied pulses to a value such that 
heating could be neglected. In order to extract the parameters determining the dynamics of the 
Josephson junction, we first used resonant activation [3T1 [TH] to extract the plasma frequency 
uo ~ 1.0 GHz at s = 0, and determine its quality factor Qo — 22. In addition, from the value 
of the plasma frequency, one can estimate the cross-over temperature between the regimes of 
thermal activation and quantum tunneling: Tco ~ fiWp/2'KkB ~ 8 mK ^ Tq. Therefore, and 
as stated above, switching does occur due to thermal activation for all experimentally accessible 
temperatures. Then, from the measurement of the number of switching events for a large number 
of current pulses, we deduced the dependence of the switching rate on the current through the 
Josephson junction. Equation ^ then yields the critical current /q ~ 0.48 yuA and the escape 
temperature Tcsc — 35 mK. The fact that the latter temperature is slightly larger than the fridge 
temperature comes mainly from the heating of the resistor Rb by the current pulses [18] . 
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Fig. 4 Switching probabilitiy Psw (left pannel) and corresponding switching rates V (right 
pannel) as a function of the reduced pulse height s = Jb/Jo for (/„i) ~ 1.1 /iA at a temperature 
To ~ 197 mK. The probabilities are calculated from statistics on 2 x 10^ to 5 x 10^ pulses in 
order to achieve a constant relative incertainty 5 x 10~*a'o on the rate F (see Appendix) , smaller 
than the size of the dots. The four curves correspond to four combinations of the signs of Vm and 
Vb- Pairs with identical sign of Vm x 14 are almost superimposed, as expected from time-reversal 
symmetry. 



For non-zero values of the average current (/„) through the NIS junction, where the statistics of 
tunneling events are expected to be Poissonian, the switching properties of the Josephson junction 
are probed as follows. Sequences of 10^ current pulses with alternating polarity and of duration 
Tp K, 50 /is are applied to the Josephson junction. In order to circumvent artifacts associated 
with possible DC offsets and thermoelectric effects, the current pulses are sent through a large, 
room temperature capacitor (200 yuF). For each pulse, the voltage across the Josephson junction 
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is monitored by two threshold detectors, each detecting the switching events in one of the two 
current directions (appearance of a positive or a negative voltage larger than a fixed threshold). 
Between each sequence of 10^ pulses, the polarity of the current (/„) through the NIS junction 
is inverted. The switching probabilities Psw can then be reconstructed for 4 combinations of the 
signs of Im and /;,. It was carefully checked that the combinations (/„ > 0, /& > 0) and {Im < 0, 
h < 0), on the one hand, and (/„ > 0, /b < 0) and (/,„ < 0, /;, > 0), on the other hand, give the 
same results, as expected from symmetry considerations. As an example, the dependence of Psw 
on s = \h\ /Iq is shown in Fig.[4]for (/„) = ±1.1 ^A. 

This measurement is then repeated for various pulse heights, and for various Im- The range 
of pulse heights was chosen to keep the relative uncertainty on the switching rates small (see 
appendix). Information on the average noise and its asymmetry are extracted from these data. 
We first focus on the average noise, which is the largest effect, seen in the escape temperature Tosc 
according to Eq. ^. 

From Ps'w, one extracts the switching rate F = ^'''p^ log(l — Psw) and the curve F(s) is fitted 
with Eq. |T]). The only fitting parameter is the escape temperature Tesc, which does not change 
for opposite current directions within the experimental error margins. The dependence of Tosc on 
{Im) is shown in Fig. El for bath temperatures Tq from 20 mK to 530 niK. 
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Fig. 5 Escape temperatures extracted using Eq. [T] (with /o = 479 nA, C = 36 pF) from the 
measurement of the rate as a function of a for several values of the current in the tunnel junction 
Im and of the refrigerator temperature Tq. The lines show the behavior expected from Eq.Q, 
in which a = e_R||/2/cs was taken as a fit parameter (aflt — 0.48K//iA). 



Apart from the 20 mK data, one finds (lines in Fig.[5|) Tcsc — 0.95 Tq + o,{Im) , with a ~ 
0.48 K/^A. This slope can be compared to the prediction ^, in which R has to be replaced with 
the parallel combination of Rf and Rb'. at ei?|| (/,„) 3> 2kBTo, one expects 



Tesc = i?||(l/i?/ + l/i?fc)To + ei?|| 
« 0.95To + ei?||(/,„)/2A:B, 



'/2kBCoth[eRt{Im)/2kBTo 



(3) 



with i?|| = {1/Rf + l/Rb + ^/Rt)^^, leading to a ~ 0.55K//iA, in reasonable agreement with 
the value observed. The slight difference in a can be attributed to an imperfect modeling of the 
impedance seen by the Josephson junction near lUp/2tt « 1 GHz, where it is the most sensitive. 
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The increased escape temperature at Tq = 20 mK is well explained by electron heating in the 
resistors Rf and Rf, in presence of the currents Im and /& [TH [22] ■ 

The switching of a Josephson junction is not only sensitive to the average noise power of current 
fluctuations, but also to higher order cumulants. The contribution of the cumulants with odd order 
yield an asymmetry of the switching rate Rr{s) = r+(.s)/r_(s) — 1, with r_|_(s) the rate when Ij, 
and Im have the same sign, r_(.s) when they have opposite signs. As discussed in the appendix, 
we acquired data only for a 0.25 < Psw < 0.987, an interval in which a good precision on Ry can 
be reached with achievable numbers N of measurement pulses. Since Tosc increases with (Im) , 
data taken at various (/,„) correspond to different ranges of s. For example, at 20 mK, s « 0.9 at 
(Im) = 0.04 ^A and s sa 0.3 at (7^) = 2 /xA. This limits the accessible range for the current /„. 

The existence of an asymmetry can be seen already in the raw data (Fig. [4|). Figure [6] shows, 
for various temperatures Tq, the dependence of the asymmetry Rr on (/„) for the two extreme 
values of rate exponents B = 9.5 and B = 11.5, corresponding to Psw ~ 0.25 and Pgw ~ 0.987, 
and to switching rates F w 6 kHz and F w 90 kHz. In each case, the two curves are very similar, 
with Ry increasing continuously with (Im) ■ 
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Fig. 6 Measured asymmetry Rr between positive and negative bias schemes as a function of 
the mean current Im through the tunnel junction. Each panel corresponds to a single refrigerator 
temperature To and each symbol corresponds to a fixed value of the rate exponent B from Eq. ([T]). 
Error bars are shown as lines (see Appendix). 

In order to compare these results with theory, we first discuss how the switching of a Josephson 
junction is affected by non-gaussian white noise. 

3 Dynamics of a Josephson junction in presence of non-Gaussian 
noise 



In this section we present a theoretical analysis of the switching of a Josephson junction out of 
its zero voltage state driven by thermal as well as weak non-Gaussian noise. According to the 
experimental situation described above (Fig. [T]), the underlying stochastic processes are assumed 
to be classical, both for the dynamics of the detector, the Josephson junction, and for the electrical 
noise originating from the mesoscopic conductor, the tunnel junction. Hence, the standard RSJ 
model applies, where now the total noise consists of the Johnson- Nyquist noise SI, in the simplest 
case produced by a resistor R in parallel to the Josephson junction, and of the mesoscopic noise 
6Im from the fluctuating part of the mesoscopic current Im = (Im) + SIm- Since the dc-component 
(7m) of the mesoscopic current can be incorporated trivially, we only focus in the sequel on its 
fluctuating part Sim- Then, according to the Kirchhoff rules we have for the currents entering the 
detector 



h + Sir 



■<" - 1 



Insunp + CjV , 



(4) 
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where <~p denotes its phase difference and the other quantities are defined according to Fig. [T] To 
proceed it is convenient to work with diniensionless quantities and to scale times with wq, energies 
with Ej. and currents with /q. This way. we arrive with (p = ^/"y^o ^"^'^ ''' ~ ^oi^ ^-t a classical 
Langevin equation for the phase dynamics, i.e., 

<^ = - sini^ - — - + Sfc + Sm(T) + sth(r) (5) 

with the quality factor Qo = RCjUq of the unbiased junction, the scaled bias current st = 
Ibl Iq, and the scaled noise sources Sth = SI/Iq and s,jj = dim/lo for thermal and non-Gaussian 
noise, respectively. These two noisy forces are assumed to be Markovian, which is an accurate 
approximation since the experiment is operated in a regime where the typical correlation times are 
much smaller than the typical time scales of the junction, e.g. the inverse of the plasma frequency. 
This situation also ensures that we can work with continuous noise processes. Hence, we have for 
the thermal noise (sth) =0 and 

(sth(rKh(0)) = — <5(r), (6) 

Wo 

where we introduced the dimensional temperature 6 = k^T^/ Ej. The non-Gaussian noise pro- 
duced by the tunnel junction is taken as Poissonian with (sm) = 0, where specifically the second 
moment reads in the high voltage/low temperature limit {eRtIm 2> 2kBTo) 

{Sn.{T)Sn,{0))^F2a^AnAr), (7) 

with Am = {lm)/Io(x, a ~ ecuo/Io and the third moment is given by 

(sm(r2)sm(Ti)s,„(0)) = F3a3A^5(r2)<5(ri) . (8) 

The prefactors F2 and F3 arc the Fano factors for the second and third cumulant. Here, for the 
tunnel junction, i^2 = -P3 = 1- It is important to realize that the average mesoscopic current 
Am can be large so that the second moment ([7]) can be comparable to or larger than that of the 
thermal noise ([6]), while the third moment is still small due to a ^ 1. In fact, as discussed above, 
the Gaussian component of the mesoscopic noise effectively adds to the thermal noise to determine 
the effective temperature of the junction 

^cff = d + 6„^ 

^ e + F,^^. (9) 

So far we have assumed that the resistance seen by the Josephson junction is exclusively given 
by the resistance R in parallel to it so that ([9]) agrees with ([2]) only in the limit Rt ^ R. In a 
description closer to the experimental situation the total resistance is RRt/{R + Rt) and we have 
to replace 0off ^^ Rt/{R + Rt) ^cff , thus arriving precisely at ([2]). Another simplification is that in 
the Langevin equation ([5]) no back-action effect from the detector to the tunnel junction has been 
taken into account. The back-action appears because when the Josephson junction switches, the 
voltage across the tunnel junction no longer coincides with the external voltage, but is reduced 
by (h/2e){Cp). Accordingly, on the one hand in the friction term in ^ Qo has to be replaced 
by Rt/{R + Rt)Qo and on the other hand the third moment ([8]) carries for finite eV/9cff an 
additional contribution from the derivative of the second moment ([7]) with respect to the voltage. 
For a detailed quantitative comparison with the experiment these modifications must be taken 
into account. 

Based on ([5]) we first give a brief account of a theoretical approach developed recently by one 
of us [23], where an analytical expression for the rate and its asymmetry due to a weak third 
cumulant of a current through a mesoscopic conductor has been derived. In a second part, direct 
numerical simulations are presented, which allow for a detailed comparison with the analytical 
predictions and may thus lay out a firm basis for the comparison between theory and experiment. 
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3.1 Analytical rate expression 

As usual in rate theory, for analytical treatments it is much more convenient to work with 
phase space probabilities rather than individual stochastic trajectories. The general problem 
for non-Gaussian noise is then that the corresponding Fokker- Planck equation (FPE) based on a 
Kramers-Moyal expansion contains diffusion coefficients up to infinite order. The basic idea for 
weak non-Gaussian noise with a leading third cumulant is thus to derive an effective, finite order 
Fokker-Planck equation. Based on a cumulant expansion of the noise generating functional such 
a generalized FPE has been derived in [23j and leads to a FPE with a momentum dependent 
diffusion term, where the momentum dependence is weighted by the third cumulant. 

The existence of a time- independent switching rate is associated with a quasi-stationary nonequi- 
librium distribution, which reduces to a thermal equilibrium in the well region and deviates from 
it only around the barrier top in the case where friction is moderate to high. This in turn means 
that the barrier height compared to the thermal energy must be sufficiently large l^U / k^T^fi S> 1 
so that switching becomes a rare event on the time scale of the plasma frequency. The typical 
thermal rate expression then looks as in ([1]) and is dominated by the exponential (activation fac- 
tor) being identical to the probability to reach the barrier top from the well bottom by a thermal 
fluctuation. Accordingly, the quasi-stationary state (the fiux state) is written as a product con- 
sisting of the thermal distribution and a 'form factor' describing the deviation from it around 
the barrier top. The latter one leads to the rate prefactor, while the former one specifies the 
activation factor. Here, we restrict ourselves to the exponential factor, which provides by far the 
dominant contribution in the present experimental set-up. Within this theoretical framework in 
[53] an analytical expression for the exponent including leading corrections due to a third cumu- 
lant have been obtained. The rate takes the form T± oc exp[— i?(l =F j^D] with the dimensionless 
barrier height B = AU/{Ej9cs) and the correction g due to the third cumulant such that r_|_ 
corresponds to Am > and r_ to A,n < 0. Note that we thus consider the bias current S(, as 
being always positive, while the situation in which one keeps A,„ > and inverts the bias current 
follows straightforwardly. 

In the dimensionless quantities introduced above one has 



Oticfi 

The rate asymmetry Rr = r+/r_ — 1 is found as 



(10) 



Rt 



exp 



16V2a3 |Am| g' (1 - Sb) 

99ls 5 + g2 yr 



5/2 



(11) 



with the quality factor of the biased junction Q = Qo(l — sl)^^^. For weak friction Q ^ 1, the 
only dependence of i?r on friction is through the effective temperature (see Eq. ^). When the 
asymmetry is small enough, it is directly proportional to the Fano factor F^, leading to a clear 
measure of the third order cumulant of the noise (Fig. Ub)- 

Let us first look at the dependence of Rr when the mesoscopic current varies (Fig. [7^) . For 
small Am the effective temperature coincides with the actual temperature ^cff ~ ^ a-nd -Rr grows 
linearly with A,„. For large Am the heating originates solely from the mesoscopic current, i.e. 
^cff ~ ^m oc Am, so that Rr decreases roughly oc 1/Am and the rate asymmetry vanishes due 
to a growing second cumulant compared to the third one. Consequently, there is a changeover 
between these two regimes with a maximum of the rate asymmetry around the range, where 
« 0m = a2AmQo/2 [sec (HI)]. 

In the experiment, the junction is operated essentially at a fixed switching rate, with Sb tuned 
at each mesoscopic current Am to a value such that B is kept constant. In this situation it is 
instructive to look for the scaling properties of the rate asymmetry as a function of the bare 
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exponent B. In general, to express RY{sb) in terms of B leads to a rather involved expression. 
Transparent formulae are found if the asymmetry is sufficiently small so that a linearization of 
the exponential in (|lip applies, for vanishing thermal noise = 0, and in the limits of strong and 
weak friction, respectively. Then, one writes for (fTTj) Rr{B) w ^„ B^ and looks for that exponent 
77 for which the prefactor ^ displays the weakest dependence on B. This way, one finds in the 
regime of weak friction that 77 sa 1.45, while for strong friction one has 77 w 5/3. 

To complete this discussion, we remark that an analytical expression for the asymmetry of the 
escape rate in the limits of low and high friction has been also derived in |25| within the framework 
of an Onsager-Machlup type of functional representation of the stochastic process combined with 
an asymptotic evaluation (high barrier, weak third cumulant) of the path integrals. The expression 
(fTT|) reduces for low Q exactly to the one specified there. For high Q both results give also precisely 
the same dependence on the physical quantities, but with slightly different numerical constants. 
Further, back-action effects according to the discussion at the end of the previous section can be 
included in pT|) in the same way as in Ref. [25] since both theoretical approaches start from the 
same Langevin equation. Hence, we conclude that, up to small corrections, both theories make 
identical predictions, where the advantage of ()lip is that it provides an analytical expression for 
all damping strengths. 

3.2 Numerical simulations 

Direct numerical simulations of the escape rate of a Josephson junction in presence of a weakly 
non-Gaussian noise are based on the Langevin equation ([5]) with the mesoscopic noise term Sm(0 
being represented as 

SmJT) = y^ aSJT - Tj) - gAm . (12) 

i 

The random variables r^ are the times at which a tunneling event occurs in the mesoscopic 
tunneling junction. Those tunneling events correspond to a Poisson process so that the intervals 
between two consecutive events St = tj+i — t^ are independently distributed according to p{Si) = 

The numerical integration of Eq.([5|) should be performed with some care as the non-Gaussian 
effects we are looking for are rather weak. We have therefore used two independent schemes and 
checked for several cases that they produced identical results within statistical errors. 

3.2.1 First scheme 

The first scheme is a first order direct integration of the equation of motion. We discretize the 
time into steps of size h. Noting ip^ = ip{T = ih) and pi = (fir = ih), we have 



with 



Pl+l - 


-P. = (f,-p,/Qo)/i + sr + Sf 






Vt+1 - 


- "ft == Pth 






F, = - 


Ai+l)h 

- sin ^i + s and T,J ~ dr sr (t) 

Jih 


(T-- 


- m; th) 



(13) 



The variables T,f^ are independent and Gaussian of variance 29/{Q[)h). In this first scheme, we 
have directly sampled the variables Si through S'i = — \og{Xi/ K^) where Xi is a random number 
uniformly distributed within [0, 1]. E™ = aN — aAm counts the number N of tunneling events in 
the window [ih, (i + l)h]. 
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Fig. 7 Numerical results versus theory, a) numerical simulations for the asymmetry of the 
tunneling rates as a function of the current Am at a fixed value of the exponent B = 6 in the 
rate ((TJ and for two values of the temperature d = 0.005 (a) and Q — 0.01 (o). Qq = 8, a — 0.02. 
The maximum occurs at the crossover Om = 0. The lines are plotted using Eq. (|ll|l with the 
same parameters, b) The lines are plotted using Eq. (|ll|l with the same parameters as in a) 
and 9 — 0.005. Each line corresponds to a different Fano factor for the third order cumulant of 
the noise. The experiment and all the other figures describe the case F3 = 1. c) asymmetry as 
a function of the exponent B for Qq = 0.5, a — 0.1, 6 — and A^ = 40. The line is plotted 
using Eq. (|lip with the same parameters, d) asymmetry as a function of the exponent B for 
Qo = 10, Q = 0.02, 9 — and Am = 40. The various symbols stand for the first integration 
scheme (o: N = 10^, h = 0.006) and second integration scheme (■: iV = 10^, h = 0.08 ;♦: 
A'^ = 10^, h = 0.0044). The line is plotted using Eq. (fTTj) with the same parameters. 



3.2.2 Second scheme 
The second scheme is a direct extension of Ref. [26] and is valid up to second order in h, 

(F, + F,+i)h/2 + S- + Sf - (^,+1 - ^J/Qo 

p,h [1 - h/{2Qo)] + F,h^/2 + (S- + Sf )/i/2 . (14) 



Pi+i - Pi 



There, we have directly sampled the number N of tunnchng events which is distributed according 

-A„,h (Am^)' 

i! 



to P{N ^i) ^e- 



3.2.3 Numerical results 

To compute the escape rate F, we integrate the equations of motions N times (typically N ~ 10^) 
and save the time Tj where (p escapes its well (taken as \(p\ > 3.5). We checked that the tj are 
distributed according to an exponential law P{tj < r) = 1 — e~^'^. Noting Tp the maximal time 
of one simulation, and < n >« e~'"'^p the measured probability of non switching, our estimator 
for F is jT = illln> ■ "^^^ asymmetry of the escape rate is defined as i?r = r(Ani)/F(— A„i) — 1. 

For the discussion of our numerical findings let us start with a comparison with the analytical 
predictions presented in the previous section. In Fig|7^ the rate asymmetry as a function of 
the average mesoscopic current Am is depicted for fixed barrier height B = 6. In accordance 
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with theory Rr shows a non-monotonous behavior with a maximum around that value of Am 
where 9 Ki 0^. For higher thermal noise the rate asymmetry shrinks for lower Am, but becomes 
independent of thermal noise in the decreasing tail for higher A,„- Apparently, while the analytical 
expression provides the overall tendency correctly, there are quantitative deviations from the 
simulated data. The reason for that will be discussed below. 

The scaling of the rate asymmetry as a function of B is shown in FiglTt for strong (Qo = 0.5) 
and in FiglTli for weaker (Qo = 10) friction. In the former case we find an accurate agreement 
with the analytical expression (jlip even for low barriers. In the latter situation the analytical 
expression describes the scaling with respect to B correctly, however, for lower barriers the absolute 
values exceed the numerical ones by about 100%, while for somewhat higher barriers, for which 
convergent data could be obtained for 6 < i? < 8, the deviations become smaller, namely 50%. 
Note, that for weaker damping the variance of the numerical data around an average behavior 
increases and is typically of the order of 10%. For even larger values of B the switching becomes 
an extremely rare event and the statistics deteriorate drastically. The obvious discrepancy for 
lower barriers that has already been seen in Figl?^ does not come as a surprise since, as known 
from standard Kramers' rate theory, escape rates are calculated as expansions in 1/B and the 
exponential activation factor dominates only for sufficiently large B. In particular in the weak 
friction regime finite barrier corrections as well as the rate prefactor tend to play an even more 
important role when one approaches the energy-diffusion limit, where the rate vanishes with 
decreasing friction strength. The turnover between energy-diffusion and spatial-diffusion occurs 
around Qo ~ ^ so that the simulated data displayed in FiglTjl in the range of lower barrier heights 
are certainly influenced by these effects. 

The overall picture emerging is hence that within the statistical uncertainties of the simulations, 
the analytical theories developed for the switching of the Josephson junction predict quantitatively 
the influence of a weak third cumulant of non-Gaussian noise up to minor deviations due to pref- 
actor effects over a broad range of parameters. As expected, larger discrepancies occur for weaker 
friction and lower barriers, a range beyond the strict applicability of the analytical treatments. 

4 Discussion 

Figure [8] presents a comparison between the measured and the predicted asymmetry in the switch- 
ing rates of the Josephson junction in presence of a current in the tunnel junction. The data of 
the left pannel (a) correspond to the intermediate the switching exponent B = 10.5. The right 
panel (b) shows the theoretical predictions Eq. (|lip for the same parameters, and assuming that 
the temperature is Tq. The disagreement is significant, both in the magnitude (asymmetries found 
in the experiment are larger than in theory) and in the current-dependence. In particular, the 
increase of the asymmetry with (/„) is in strong contradiction with theory, which predicts a max- 
imum when the shot noise due to the tunnel junction is similar in magnitude to the thermal noise 
seen at zero current. Two reasons can be invoked for this disagreement: an imperfect modeling 
of the experiment by theory, or experimental artifacts. 

4.1 Effects not considered by theory 

The experimental setup presents differences from the simplified circuit described by theory: the 
RC combination that allows in the experiment to separate low and high frequency fluctuations is 
not treated by theory; the resistance Rf in series with the tunnel junction also produces noise that 
modifies the voltage across the tunnel junction, a configuration in which feed-back effects lead to 
an asymmetry of the current fluctuations [8] that has not been considered; the voltage across the 
Josephson junction has fluctuations due to plasma oscillations, leading to voltage fluctuations that 
have been neglected. At least this last point can be discarded, since the corresponding voltages 
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Fig. 8 a) Measured asymmetry Rr as function of the current /,„ at a fixed value of the 
exponent B = 10.5 in the rate ^ for various temperature, b) Same quantity calculated using 
(|lip . assuming the actual temperature is the refrigerator temperature. 



Si ifiQUJp arc in the fiV range. More work in necessary to analyze the effects of the two other main 
differences between experiment and theoretical modeling. 
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Fig. 9 a) Plot of the asymmetry Rr for a fixed value B — 10.5 assuming that the actual value 
of the current s is given by s <— s ± \{Im)/Io because of a possible leak through the capacitor 
Ci. Besides, the temperature 6cS was extracted from Fig[5] instead of being calculated from 
Eq.®. Here, A — 5.5 x 10~* leads to the best fit with the data, b) Same thing assuming that 
the third moment is 0, that is <; = 0. Here, A = 7 x lO"** leads to the best fit with the data. 
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4.2 Experimental artifact? 

We have also explored another possible source of asymmetry associated with an hypothetic im- 
perfection in the experiment, namely a leak across one of the coupling capacitances C or Ci or 
a touch to ground of the corresponding bias line. In this case, a fraction A of the DC current 
through the junction passes through the Josephson junction instead of returning entirely through 
the resistor Rf-. at {Im) > 0, the switching rate at /& > is increased, whereas at at (/„) < it 
is decreased. The corresponding asymmetry can be evaluated using Eq. ([T]) with B evaluated at 
s = (/(, ± A {!,„)) /Iq- For A <C 1, one obtains 

Rr~3B{l-.s)-'X{Ira)/Io- (15) 

The best fit of the data with this expression in which the measured escape temperatures are 
plugged in, obtained with A = 7 x 10"'^, is presented in Fig. [Ob. The general trend is well 
reproduced, but systematic deviations are observed for the data at Tq = 20 niK and at currents 
above 1 fiA. In Fig. [9^. we fitted the data with the conjugated effect of a leak and of the 
asymmetry expected from the third order cumulant of the current fluctuations (replacing again 
the effective temperature ^cff by the measured escape temperature). In the best fit, obtained with 
A = 5.5 X 10""', all the data are precisely accounted for. If this scenario is the correct one, the 
asymmetry due to the tunnel junction is clearly seen in the data at the lowest bath temperature, 
but the uncertainty is too large for other temperatures. The value of A corresponds to a leak 
resistance of Rt/X « 3 MH. {Rf/X « 0.4 Mfl) if the leak comes from the capacitor Ci (C). The 
capacitors were tested after fabrication, but a leak might have appeared afterwards. A touch to 
ground of the corresponding connection line at this intermediate value seems unlikely. 

5 Conclusions 

In conclusion, we have measured the switching rate of a Josephson junction when coupled to a 
voltage-biased NIS junction. The main effect is a well explained increase of the escape temperature 
due to shot noise through the junction. The asymmetry of the rates for opposite signs of the bias 
of the NIS junction is larger than expected from the dynamics of the Josephson junction in 
presence of non-Gaussian noise. Part of the discrepancy might be due to an oversimplification in 
the modelling of the circuit. It is however suspected that a current leak through a capacitor, is 
responsible for a large fraction of the signal. Experiments less prone to experimental artifacts, 
with an improved decoupling of the detector from the noise source at DC, are in progress. 

A Appendix: Sensitivity of the switching probabihty to a rate 
asymmetry 

For each current pulse we recorded the boolean 1 or depending whether the voltage across 
the Josephson junction crossed the threshold or not. The booleans arc N independent random 
variables Xi with binomial distribution, Psw being the expectation for Xi = 1. Therefore, the 
standard deviation of the average N"^ J2i=i -^i i^ \/ Psw{^ ~ Psw)/N. And using 

P,^ = l-e-^"p, (16) 

one gets the standard deviation of the switching rate 



''^^w^^r-'- <"> 
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Therefore, the uncertainty on the measurement of i?r is given by (see Eq. ([T7|)): 



^""^ -^T~ -^im^^)W^- ^''^ 



This function presents a minimum at dRr ~ 2.5/\/N for Psw ~ 0.8, and grows rapidly near 
Psw = and Psw = 1- In the experiment, we chose the range of s so that 0.25 < Psw < 0.987 (or 
equivalently 9.3 < B < 12 in the notations of Eq. (IT|))for which SRr < 4/\/]V', and adapted N at 
each value of Psw to keep SRr at a fixed value. 
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